VARIABLE STRUCTURE ALGORITHM APPLYING 3-STAGE METHODS                                      OF RUNGE-KUTTA AND ROSENBROCK TYPES

Novikov Evgeniy Aleksandrovich, Doctor of physical and mathematical sciences, professor, chief researcher, Institute of computational modeling of the Siberian branch of RAS (building 44, 50 Akademgorodok street, Krasnoyarsk, Russia), atp@sstu.ru

519.622

Background. The Cauchy problem for large-scale stiff systems arises in simula-tion of physical and chemical processes, in approximation of partial differential equations by a system of ordinary differential equations and in plenty of other impor-tant applications. Taking into consideration a large number of factors in model de-velopment leads to expansion of a class of problems, determined by stiff systems of high dimension. The complexity of practical problems leads to additional require-ments to computational algorithms.
Materials and methods. Dealing with high dimensionality of a stiff system of or-dinary differential equations, the main computational expenses concern the Jacobi matrix decomposition. In some algorithms one may use freezing of the Jacobi matrix, i.e. applying the same matrix over several integration steps. The problem of freezing is solved rather easy in those methods, the stages of which are computed using some iterative processes. For non-iterative numerical formulas the freezing is quite a diffi-cult problem. In this study, the costs reduction was achieved by combining explicit and L-stable methods using the criterion of stability in calculations.
Results. The author has created an algorithm of variable structure integration, based on the explicit scheme of the Runge-Kutta type and the L-stable method of the Rosenbrock type. Both schemes have the third order of accuracy. An efficient nu-merical formula was chosen according to the criterion of stability at each step of in-tegration. Estimation of maximum eigen value, which is necessary to switch between the methods, for explicit numerical schemes was determined by power iterations us-ing already computed stages, and using the Jacobi matrix norm for the Rosenbrock type method. The researcher also formulated inequalities for accuracy and stability control. The article adduces the results of calculations.
Conclusions. The integration algorithm is aimed at solving stiff problems of high dimension. Numerical results confirm the efficiency of the constructed algorithm.

stiff system, schemes of Runge-Kutta and Rosenbrock types, accuracy and stability control, automatic selection of a method.

1. Hairer E., Norsett S. P., Wanner G. Solving Ordinary Differential Equations I: Nonstiff and Differential-Algebraic Problems. Berlin: Springer-Verlag, 1993, 528 p.
2. Hairer E., Wanner G. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Berlin: Springer-Verlag, 1996, 614 p.
3. Novikov E. A., Shornikov Yu. V. Komp'yuternoe modelirovanie zhestkikh gibridnykh sistem [Computer simulation of stiff hybrid systems]. Novosibirsk: Izd-vo NGTU, 2012, 451 p.
4. Rosenbrock H. H. Computer. 1963, no. 5, pp. 329–330.
5. Novikov E. A., Novikov V. A., Yumatova L. A. Zhurnal vychislitel'noy matematikt i matematicheskoy fikhiki [Journal of calculus mathematics and mathematical physics]. 1987, vol. 27, no. 3, pp. 385–390.
6. Novikov E. A., Novikov A. E. Mathematical Models and Computer Simulations. 2010, vol. 2,no.4,pp.443–452.
7. Novikov E. A. Soviet reports on mathematics. 1984, vol. 30, no. 2, pp. 358–361.
8. Novikov E. A. Yavnye metody dlya zhestkikh sistem [Explicit methods for stiff systems]. Novosibirsk: Nauka, 1997, 198 p.
9. Fehlberg E. Low order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems. Washington: NASA T.R.R., 1969, 316 p.
10. Dahlquist G. A special stability problem for linear multistep methods. BIT. 1963, vol. 3, pp. 23–43.
11. Demidov G. V., Yumatova L. A. Issledovanie tochnosti neyavnykh odnoshagovykh metodov [Accuracy research for non-explicit one-step methods]. Novosibirsk: Izd-vo VTs SO AN SSSR, 1976, (Preprint № 25).
12. Demidov G. V., Novikov E. A. Chislennye metody mekhaniki sploshnoy sredy [Numerical methods of continuum mechanics]. 1985, vol. 16, no. 1, pp. 27–39.
13. Enright W. H., Hull E. Comparing numerical methods for the solutions of systems of ODE’s. BIT. 1975, vol. 15, pp. 10–48.